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ABSTRACT 


This  investigation  examines  the  inherent  potential  flows 
and  laminar  boundary  layer  characteristics  on  the  nose  surface 
of  axisymmetric  bodies  that  can  postpone  the  natural  transi¬ 
tion  process  from  a  laminar  to  turbulent  flow.  Based  on  this 
analysis,  an  effort  is  made  to  find  a  correlation  between  the 
geometrical  features  of  a  given  body  shape,  and  the  ability  of 
that  shape  to  delay  transition.  Ellipsoidal  forebodies  for 
which  stability  computations  have  been  performed,  are  pre¬ 
sented  to  show  those  attributes  of  a  good  nose  design.  These 
examples  give  some  direction  to  future  work  in  obtaining  the 
optimal  design. 

Prior  to  the  discussion  of  the  effects  of  nose  shape  on 
transition,  a  brief  review  of  the  method  commonly  used  for 
predicting  transition  is  presented  followed  by  a  discussion  of 
the  distinctive  features  of  axisymmetric  flows  as  compared  to 
the  widely  discussed  two-dimensional  flows.  This  is  essential 
background  material  and  forms  the  basis  of  most  of  the  argu¬ 
ments  presented  on  the  choice  of  a  good  bow  shape  design. 


ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  is  funded  by  the  Naval  Sea  Systems  Command 
(SEA  037)  under  the  FDT&E  Exploratory  Development  Program. 


1.0  METHOD  FOR  PREDICTING  TRANSITION 


1.1  BACKGROUND 

The  development  of  the  viscous  boundary  layer  on  the  nose  of  an  axi- 
symmetric  body  normally  starts  at  the  stagnation  point  as  a  steady  lami¬ 
nar  flow  growing  in  thickness  along  the  downstream  direction  until  irregu¬ 
lar  fluctuations  or  bursts  appear  causing  the  process  of  transition  from 
the  laminar  to  turbulent  flow^".  The  transition  process  is  caused  by  "in¬ 
stability"  of  the  flow,  that  is,  the  tendency  for  small  natural  distur¬ 
bances  (due  to  noise,  mechanical  vibrations,  surface  roughness,  non-uni¬ 
formity  of  the  oncoming  stream,  etc.)  to  be  amplified  into  substantial 
fluctuations . 

The  classical  method  for  predicting  the  location  of  transition  is 
based  on  stability  theory  which  attempts  to  predict  if  these  natural 
disturbances,  inherent  in  all  flows,  will  grow  or  decay  during  the  down¬ 
stream  development  of  the  laminar  boundary  layer. 

The  analysis  assumes  that  these  disturbances  are  represented  as 
infinitesimal  two-dimensional  velocity  fluctuations 

u'  "  <P(y)  exp  [i (ax-cot)] 


W  J 


O  * 

ms 


where  a  =  +  ia^  -  complex  wave  number  of  disturbance 

c  =  2Tt(o/ar,  phase  velocity  of  wave 
co,  circular  frequency  of  wave 

Instability  of  the  flow  occurs  when  the  spatial  amplification  rate 
is  positive  or  <  0,  for  which  the  disturbance  begins  to  grow  spatially 
in  a  complex  three-dimensional  manner  until  the  eruption  of  a  turbulent 


VW>. 


burst.  The  growth  pattern  of  these  discrete  unstable  waves  can  be 


LDeMetz,  F.C.  and  M.J.  Casarella  "An  Experimental  Study  of  the  Intermit¬ 
tent  Properties  of  the  Boundary  Layer  Pressure  Field  During  Transition 
on  a  Flat  Plate,"  NSRDC  Report  4140  (November  1973). 

2 


characterized  by  the  spatial  amplification  ratio,  ln(yy  ),  which  measures 

o 

the  amplification  from  the  initial  instability  (neutral  stability  point) 
to  the  onset  of  burst  eruptions.  The  Smith-Gamberoni  method  of  predicting 

transition  is  based  on  the  criterion  that  the  onset  of  bursting  will  occur 

A  2 

when  ln(—  )  :  9-11.  However,  the  pioneering  work  by  Granville  was  the 
n 

o 

first  attempt  at  correlating  the  conditions  of  the  flow  instability  with 
the  onset  of  fully  turbulent  flow.  His  work  identified  the  inherent  flow 
characteristic  that  allow  realistic  predictions  for  a  wide  class  of  body 
shapes. 

It  has  been  shown  experimentally  that  these  two-dimensional  regular 
waves  do  not  properly  represent  the  growth  of  the  waves  on  axisymmetric 
bodies  prior  to  the  breakdown  of  the  flow  into  turbulent  burst.  Further¬ 
more,  the  amplification  ratio  concept  of  predicting  transition  is  based 
strictly  on  linear  stability  theory  for  infinitesimal  disturbances  and 
does  not  account  for  the  actual  phenomena  of  three-dimensionality,  non¬ 
linear  interaction  of  the  waves,  and  vortex  breakdown.  However,  the 
Smith-Gamberoni  method  appears  to  predict  transition  within  a  20  per  cent 
range  for  low  Reynolds  number. 

In  order  to  predict  the  growth  of  the  wave  by  this  method  and  its 
subsequent  amplification  rate  when  ct^  <  0,  numerical  solutions  of  the 
eigenvalue  problem  posed  by  the  Orr-Sommerf eld  equations  are  required  at 
each  station  along  the  axisymmetric  body  surface.  This  requires  a 
significant  amount  of  computational  effort.  The  input  for  this  program 
requires  knowing  the  laminar  boundary  layer  velocity  profile  (yy  )  and 


^Granville,  P.S.,  "The  Calculation  of  the  Viscous  Drag  of  Bodies  of 
Revolution,"  DTMB  Report  No.  849  (July  1953). 
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the  local  displacement  Reynolds  number  R. .  at  each  station  along  the 

6  3 

axisymmetr ic  body  surface  prior  to  the  onset  of  transition. 

The  Smith-Gamberoni  method  of  prediction,  based  on  solutions  of 
the  Orr-Sommerfeld  equation,  is  widely  used  because  it  reflects  the  best 
current  state-of-the-art  technique.  The  results  appear  to  have  some 
correlation  with  experimental  results  obtained  in  model  testing  at 
Reynolds  numbers  (R^  =  8-10x10^).  These  values  are  significantly  lower 
than  those  existing  in  an  operational  environment  (e.g.,  R^  *  50  - 
110  x  10^).  Experimental  tests  at  higher  Reynolds  numbers,  which  are 
lacking,  are  required  in  order  to  gain  confidence  in  the  validity  of  the 
prediction  technique  for  even  the  most  conventional  shaped  bodies.  How¬ 
ever,  the  Smith-Gamberoni  method  can  be  used  to  make  comparative  evalua¬ 
tions  of  various  shapes  in  order  to  gain  insight  into  the  so-called  "best 
preliminary  design."  This  is  essentially  the  philosophy  underlying  the 
remarks  contained  in  this  report. 

1.2  SUMMARY  OF  COMPUTATIONAL  PROCEDURE 

The  overall  algorithm  for  predicting  the  transition  point  is  shown 
in  Figure  1.  Four  distinct  phases  of  the  computational  algorithm  are 
required.  These  are  shown  schematically  in  Figure  1  and  will  be  briefly 
discussed.  Appendix  A,  which  is  a  condensed  version  of  the  equations 
derived  in  Reference  3,  contains  the  derivation  and  summary  of  the 
equation  used  in  this  discussion. 

A .  Axisymmet r ic  Potential  Flow  Computation 

U 

e 

By  using  the  Neumann  Program,  the  pressure  distribution  C  and  — 


Niedzwecki,  J.,  "Laminar  and  Turbulent  Incompressible  Boundary  Layers  on 
Bodies  of  Revolution  in  Axial  Flow,"  Doctoral  Dissertation,  Catholic 
University  of  America,  Dept,  of  Mechanical  Engineering  (Feb  1977). 
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are  computed.  From  these  results  one  must  numerically  evaluate  M,  given 


by  equation  (A5)  as 


where 


2s3  dU 

M3(s3) - — - 

Ue  ds3 


'1  +  (r^)  dx 


and  s3  is  the  distance  along  the  body  surface,  rQ(x)  the  radius  of  the 

body,  and  x  is  the  axial  distance  measured  from  the  nose  of  the  body. 

U 

8  m 

In  the  vicinity  of  the  stagnation  point  of  a  blunt  body  —  «  s3  ,  hence 

oo 

lim  (M3)  ■  2m.  For  the  axisymmetric  stagnation  flow  on  a  disk,  m  ■  1, 
s3-*0 

hence  lim  M3  =  2.0.  This  value  is  used  in  all  the  computer  studies. 
s3**0 

In  general  2.0  _>  M3  _>  0  and  decreases  with  increasing  s3  in  the  favorable 
pressure  gradient  region. 

B.  Maneler  Transformation  into  Two-Dimensional  Coordinates 


The  axisymmetric  coordinates  (s3,  y3)  and  the  pressure  gradient  para¬ 
meter  M3  are  transformed  into  the  equivalent  two-dimensional  similarity 
coordinates  (C,n)*  and  pressure  gradient  parameter  B  through  the  equations 

(A2),  (A3),  and  (A4)  given  as 

2 _  _ 2  s  3 

C  =  u  U  Rp,r  9T3 
e  e  D  o  2P 

- 

n  =  y  — — Y3 

”  s3 


M3 

6  f,  =  - ? - 

2P3(s3) 


where  s3  =  s3/D,  y3  =  y3/D,  r  =  r  /D,  U  -  U  /Ut 


*These  are  called  the  Mangler-Levy-Lees  transformations,  see  Appendix  A. 


2  JU  r^ds 

J0  6  °  3  ■ 

The  local  value  of  6  along  the  nose  surface  is  the  single  most  important 
external  flow  characteristic  governing  the  laminar  boundary  layer  pro- 
perties  and  consequently,  affecting  the  location  of  transition.  There¬ 
fore,  some  comments  are  in  order  regarding  the  magnitude  of  this  function 

in  the  favorable  pressure  gradient  region  downstream  of  the  nose  stagnation 

U 

0  m  p 

point.  Again,  assume  —  “  s  and  r  «  sr  for  the  stagnation  region  on 

U  3  o  3 

CO  J 

the  axiymmetric  body,  then  lim  (M  )  =  2m  and  lim  (2P2)  ■  2p+m+l  hence, 

s,-*-  0  3  s  -K)  3 

3  3 

6(0)  =  »  - —  .  For  the  disk  stagnation  flow,  m  -  1  and  p  =  1,  there- 

Zp-rnrrl 

fore  6(0)  =  0.5  and  this  value  is  used  as  the  starting  value  for  the 
computer  solutions.  In  general,  the  range  of  values  in  the  favorable 
pressure  gradient  region  downstream  of  the  nose  are 


2  i 

2.0  >  M  (s  )  >  0  2  >  p  (s  )  >  - 

“  3  3  -  “  3  3  /2 


0.5  >  6(0  >  0. 


C.  Two-Dimensional  Laminar  Boundary  Laver  Computations 


The  two-dimensional  laminar  boundary  equation  in  the  transformed 
4 

coordinates  is  given  as  , 


f'"  +  f"  +  8[  1-  (f '  )  2  ]  -  2£  [  f '  "  f"  If] 

at,  dt. 


*It  has  been  brought  to  our  attention  by  Dr.  von  Kerczex  that  B>1  can  occur 
on  flat  faced  bodies  of  revolution  downstream  of  the  stagnation  point. 


where 


—  =  f'(n)  and  8  =  function  (M  ,P  ,£)• 

Ue  3  3 

This  equation  is  numerically  solved  for  the  velocity  profile  ^  using  the 

4  e 

Cebeci-Smith  method  .  The  displacement  thickness  6*  along  the  body 

3 

surface  is  then  computed  through  the  equations  (A8  )  and  (A9  ) 


where 


s  U 
3  e 


P3 


63 


U  6* 
l  =— e_3- 

<5*  v 


U  D 


and 


r 


<5*  (6,0  =  J  (1-f * )  dn. 

o 


From  this  equation,  it  is  seen  that  the  body  Reynolds  number  Rp  has 

a  direct  effect  on  the  displacement  thickness  5*  along  with  the  value  of 

3 

6  which  enters  through  the  computations  of  6*(B,0- 
D.  Orr-Sommerfeld  Stability  Computations 

Using  the  velocity  profile  77  and  the  displacement  thickness  6*  in 

u  -3 

e  J 

the  axisymmetric  coordinates  (s  ,y  ),  the  Orr-Sommerfeld  equation  is 

3  3 

solved  numerically  at  each  point  along  the  body  surface  s  for  selected 

3 

frequencies  of  the  disturbance  waves.  When  one  or  more  of  these  dis¬ 
turbance  waves  are  found  to  grow  (a  <0),  then  the  amplification  rate 


8 


a 


ln(x  >  is  computed  for  each  of  these  growing  waves  until  at  least  one 


wave  is  found  to  supposedly  erupt  into  a  turbulent  burst  (In  —  :  9-11) . 

o 


This  location  is  then  defined  as  the  location  of  transition  St  where  St 


funct  (shape,  R^) . 

The  solution  of  this  equation  is  a  computational  task  of  order  of 
magnitude  greater  than  the  three  other  tasks.  A  major  research  effort 


is  being  pursued  in  developing  more  efficient  algorithms  that  will  reduce 

^  . _ _  .  t  .  I  .  1  I  TT  .  r  l  5  ,  1  J 


the  computer  execution  time.  The  recent  work  by  Von  Kerczek  indicates 
that  this  task  can  be  performed  within  reasonable  execution  times  and  cost. 

Extensive  numerical  solutions  to  these  equations  have  been  obtained 
for  a  wide  class  of  laminar  velocity  profiles  referred  to  as  similar  flows. 


These  are  the  cases  when  8  =  constant  and  f  =  f(n)  and  includes  conical 

1  6  .  I  .  .  .  1  i  fT<1  T.  T  T  1  .1  . 


shape  noses  of  half  angle  <p .  These  results  clearly  show  the  strong 
dependence  of  the  location  of  instability  on  the  values  of  8.  More 


specifically,  these  results,  to  a  rough  approximation,  give 


Similar  Flows 


a 

* 

8 

(R5*)crit 

6 

<j>(cone  half  angle) 

ohm 

1.0 

12,490 

.500 

90. C° 

5 

0.8 

10,920 

.421 

81.6° 

MM 

0.6 

8,890 

.333 

71.3° 

is 

5 

0.4 

6,230 

.210 

58.2° 

0.2 

2,830 

.125 

40.3° 

IK 

0.0 

520 

.064 

27.7° 

.033 

19.0° 

Von  Kerczek,  C.,  Private  Communication. 


Evans,  H.L.,  "Laminar  Boundary-Layer  Theory,"  Addison-Wesley  Publishing 
Co.,  London,  p,121,p.70  (1968). 


White,  F.M.,  "Viscous  Fluid  Flow,"  McGraw-Hill,  New  York,  p.  341,  p.  405 
(1974). 
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2.0  DISTINCT  FEATURES  OF  AXISYMMETRIC  FLOWS 


2.1  DISTINCTION  BETWEEN  TWO-DIMENSIONAL  AND  AXISYMMETRIC  FLOWS 


A  clear  distinction  should  be  made  between  the  axlsymmetric  coordinates 

(s  ,y  )  where  s  is  the  distance  along  the  nose  surface  and  y  is  normal 

3  3  3  3 

to  the  surface,  and  an  equivalent  two-dimensional  system  (s  ,y  )  such  as 

2  2 

that  along  a  flat  plate.  Figure  2  summarizes  the  transformations  and 
boundary  layer  characteristics  for  each  of  these  cases. 

The  two  systems  are  considered  equivalent  if,  after  the  transformation, 

the  transformed  variables  £,n  are  equivalent  and  6  ■  B  .  In  order  to 

2  3 

make  a  comparison  between  the  two-dimensional  and  axlsymmetric  system, 
two  cases  can  be  equated: 

CASE  I  -  Equivalent  Downstream  Distance 

Given:  U  (s  )  and  r  (s  ) 

e3  3  °  3 

Assume:  s  =  s 

3  2 

Then,  the  two  flows  are  equivalent  when 


y 

2 


r  y 
°  3 


r2  U  (s  )  . 
o  e  3 


CASE  II  -  Equivalent  Pressure  Gradient  in  External  Flow 

Given:  U  (s  )  and  r  (s  ) 

e3  3  0  3 

Assume:  U  (s  )  =  U  (s  )  or  M  *  M 

e2  2  e3  3  3  2 

Then,  the  two  flows  are  equivalent  when 
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Figure  2  -  Comparison  of  Two-Dimensional  and  Axisymmetric 
Boundary-Layer  Properties 


It  is  erroneous  to  make  a  direct  comparison  of  the  axisymmetric  flow 


field  with  those  of  a  two-dimensional  system  in  which  one  equates  both 
the  distance  downstream  along  a  flat  plate  with  the  axisymmetric  surface 

distance  (s  =  s  )  and  assumes  identical  external  pressure  gradients 

2  3 

(M  =  M  ).  The  values  of  8  used  for  the  computation  in  the  transformed 

2  3 

coordinates  (£,n)  will  be  different  and  thus  the  computed  boundary  layer 
characteristics  are  not  identical.  By  examining  the  equations  given  in 


Figure  2,  one  can  show  that  in  a  favorable  pressure  gradient  region  near 
the  nose  or  leading  edge 


two-dimensional  flow 


axisymmetric  flow 


1.0  >6  >0 

2 


0.5  >  8  >  0 

_  2  - 


Consequently,  the  pressure  gradient  parameter  8  on  the  equivalent  flat 

2 

plate  will  be  higher  by  a  factor  of  2  from  the  true  axisymmetric  pressure 
gradient  6^  despite  the  fact  that  both  systems  have  the  same  external 
flow  fields. 

It  can  also  be  shown  that  the  displacement  thickness  ratio  between 

the  actual  6*  and  the  equivalent  two-dimensional  6*  will  be 

3  2 


6*  p  (s  )  6*(B  ,£  )  .  6*(B  ,5  ) 

_1  =  _2 2_  2  ■  2  =  1  ..  2 2_ 

6*  W  6*  (8  ,£  )  ST  «*  (6  ,£  ) 

2  3  3  3  3 


where  8  >8  and  £  >  £  . 

2  3  2  3 

Similarly,  solutions  of  the  boundary  layer  equations,  <5*(8), 
decrease  slightly  with  increasing  values  of  8,  as  shown  below'7, 


AJJ 


Elliptical  Forebodies 


t-1. 

,00 

1  t=. 500 

t= 

Similarity 

6*(g) 

<5*(8,0 

<S*(8,0 

<5*(e,£) 
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2.00 

.804 

2.00 

.804 
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.836 

1.73 

.842 
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- 

.911 

.871 

1.58 

.889 

1.66 

- 

.984 

.911 

1.47 

.949 

1.46 

.945 
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.95  8 

1.39 

1.020 

1.23 

1.040 

1.217 

.996 

1.33 

1.113 

1.06 

1.172 

-1.54 

-1.30 


"t"  in  the  above  set  of  values,  is  the  ratio  of  the  minor-to-major  axis 

5*(e2> 


Therefore 


6*(62)  <  <$*(63)  or 


0  .82. 


6*(B3) 

By  examining  various  numerical  examples,  one  obtains  the  approximation 
that 


6*  *  0.86  6$ 

3 


and  leads  one  to  speculate  that  in  general,  the  displacement  thickness 
along  an  axisymmetric  body  would  be  slightly  thinner  than  the  equivalent 
displacement  thickness  along  a  flat  plate  with  the  identical  pressure 
distribution. 

In  summary,  two  generalized  remarks  can  be  made*: 
i)  The  pressure  gradient  imposed  by  the  external  flow  on  an 
axisymmetric  body  is  weaker  than  the  equivalent  pressure  gradient  imposed 
on  a  two-dimensional  flat  plate. 


^Similar  conclusions  have  been  expressed  by  Kaups  on  page  10  of  Douglas 
Report  No.  MDC  J6530  (April  1974). 


.v  *-  V/7V 
■-  -'.••V-V.-  V.vj  - 


•  • 

.'•v- 


ii)  The  displacement  thickness  on  the  axlsymmetrlc  body,  6*  ,  is 
thinner  than  that  which  would  exist  on  the  flat  plate  with  an  equivalent 
pressure  gradient. 

The  fact  that  the  boundary  layer  characteristics  on  an  axisymmetric 
body  are  different  than  those  on  a  flat  plate  requires  one  to  interpret 
this  distinction  from  the  standpoint  of  their  effect  on  transition.  By 
comparing  the  location  of  transition  on  an  axisymmetric  body  with  that 
of  a  flat  plat  with  identical  pressure  distribution,  the  weaker  pres¬ 
sure  gradient,  in  terms  of  8,  is  an  adverse  effect  to  enhancing  a  stable 

g 

flow  while  the  thinner  boundary  layer  is  a  favorable  effect  . 

2.2  LOCALIZED  FLOW  PROPERTIES  GOVERNING  TRANSITION 

In  the  discussion  on  the  stability  prediction  method,  it  has  been 
emphasized  that  the  local  values  of  the  pressure  gradient  parameter  8 
have  a  strong  influence  on  the  location  of  transition.  Specif ically , 

that  transition  is  delayed  by  a  strong  favorable  pressure  gradient  (8 
as  large  as  possible).  The  values  of  8  are  inherently  related  to  the 
geometrical  shape  of  the  body  and,  as  previously  discussed,  require 
different  interpretations  for  two-dimensional  and  axisymmetric  flows. 

Let  us  try  to  see,  more  specifically,  how  the  magnitude  of  8 
affects  the  onset  of  transition.  Consider  a  favorable  pressure  gradient 
0.5  _>  B(s3)  ■>  0  representing  the  region  downstream  of  the  stagnation 
point.  One  can  present  a  somewhat  simplified  and  qualitative  interpreta¬ 
tion  of  the  transition  process  by  arguing  that  the  instability  of  a  dis¬ 
turbance  wave  will  occur  when  the  local  displacement  thickness  R^*(s3), 

Q 

Kaups,  K.,  "Transition  Prediction  on  Bodies  of  Revolution,"  Report 
No.  MDC  J6530,  Douglas  Aircraft  Co.,  Long  Beach,  California  (Apr  1974). 


computed  at  each  downstream  position,  exceeds  a  certain  critical  value 


(R^^crit*  Some  justification  for  using  the  dimensionless  displacement 

thickness  as  a  measure  of  the  stability  threshold  are  presented  by 
9 

Lighthill  in  which  R^  is  interpreted  as  a  ratio  of  the  convection  rate 
to  diffusion  rate  of  a  disturbance  in  the  vorticity  field.  Alternatively, 
one  could  use  the  dimensionless  momentum  thickness  R.  since  this  appears 
to  be  the  fundamental  characteristic  proposed  by  Granville  in  Reference  2. 


/u  6*  \ 

The  displacement  thickness  R^  *  I  e  3  I  increases  along  the  down- 

3  \  v  / 

stream  distance  because  of  the  growth  of  the  boundary  layer.  The  particu¬ 
lar  rate  of  the  growth  depends  on  the  local  values  of  the  pressure  gradient 
parameter  6,  and  the  body  Reynolds  number  R^.  On  the  other  hand,  the 
critical  Reynolds  number  has  an  extremely  large  value  at  the 

stagnation  point  but  decreases,  rather  drastically,  along  the  downstream 

distance  because  of  the  decreasing  8.  When  the  increasing  value  of  R^ 

3 

equals  the  decreasing  critical  value  (^<5*)crit>  transition  will  occur. 

Figure  3  schematically  illustrates  this  condition  for  three  arbitrary 

body  shapes.  The  key  to  delaying  transition  is  then  two-fold: 

i)  Extended  Favorable  Pressure  Gradient  Region 

The  local  pressure  gradient  parameter  B  should  be  as  large 

as  possible  over  an  extended  region  in  order  for  a  high  critical  Reynolds 

number  or  (R..)  =  funct  (B)  should  be  large. 

o*  crlt 

ii)  Thin  Boundary  Layer  During  Downstream  Development 

The  boundary  layer  growth  should  be  as  thin  as  possible 

during  its  downstream  development,  or  R^*  =  function  (P^BjR^)  should  be 
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Q 
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England,  p.  88-90  (1963). 
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Figure  3  -  Schematic  Diagrams  Illustrating  the  Effects  of 
Nose  Shapes  on  the  Location  of  Transition 


small.  Note  that  equation  (A9)  gives 


r  1 


-3^ 


where 


<S*(.5,0) 

vy 


.804 

1.414 
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Both  of  these  aspects  are  enhanced  by  a  B (S ^ )  function  which  decreases 
at  the  slowest  possible  rate  from  its  maximum  value  at  the  stagnation  point 

It  is  important  to  note  that  the  rate  of  growth  of  the  boundary  layer 

U  D 
00 

thickness  is  strongly  influenced  by  the  body  Reynolds  number  R^  ■  -  , 

v 

more  so  than  by  the  magnitude  of  6.  Therefore,  despite  the  advantages  of 
the  favorable  pressure  gradient  in  delaying  transition  past  the  (Cp)m^n 
location,  transition  can  still  move  forward  and  occur  within  the  favorable 
pressure  gradient  region  at  very  high  Reynolds  number  R^.  This  is  illu¬ 
strated  in  Figure  4  for  body  B  of  previous  figure. 

A  comparative  evaluation  of  several  body  shapes  can  be  made  by 
examining  the  behavior  of  the  6  curves  downstream  of  the  stagnation  point 
for  each  of  these  shapes.  These  results  can  be  obtained  from  the  potential 
flow  computation  (Program  I).  Figure  5  illustrates  these  8  computations 
for  the  standard  ellipsoidal  forebodies  obtained  from  the  Neumann  Program. 
Imposed  on  these  figures  is  the  8  curve  for  an  "idealized"  shape  body 
possessing  the  property  of  a  large  value  of  8  and  consequently  a  large 
value  of  (R,5*)cr^t  downstream  of  the  stagnation  point. 

The  growth  rate  of  the  boundary  layer  downstream  of  the  stagnation 
point  is  also  an  important  feature  in  evaluating  several  body  shapes.  The 
laminar  boundary  layer  computations  for  each  of  the  body  shapes  are 
required  (Program  II).  However,  as  one  would  expect,  the  shape  of  these 


curves  downstream  of  the  stagnation  point  are  Initially  quite  similar 


Figure  4  -  Schematic  Diagram  Illustrationg  the  Effect  of  Body 
Reynolds  Number  on  the  Location  of  Transition 


until  3  -*•  0  for  each  of  the  respective  shapes.  Figure  6  illustrates 


these  computations  for  ellipsoidal  nose  shapes  at  fixed 
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3.0  EFFECTS  OF  GEOMETRIC  SHAPE  ON  TRANSITION 
3.1  COMPARISON  OF  VARIOUS  SHAPES  AT  FIXED  R^ 


As  previously  stated,  the  first  task  in  any  evaluation  of  the 
location  of  transition  on  an  axisymmetric  nose  is  to  compute  the  potential 
flow  or  pressure  distribution  around  the  body  surface.  The  relationship 


between  the  geometrical  shape  of  an  axisymmetric  body  and  the  shape  of 


the  pressure  distribution  curve  (C  versus  — )  or  the  freestream  velocity 
Ue  X  P 

curve  (—  versus  — )  has  been  extensively  studied.  The  two  most  char¬ 


acteristic  features  are  the  magnitude  and  location  of  (C  )  ,  or  alternately 

u  uep  min 

the  magnitude  and  location  of  (~)  where  C  =  1  -  (—  )2.  Figure  7 

U  max  p  U 

U  CO  r  oo 

g 

illustrates  the  —  curves  for  a  family  of  elliptic  noses.  The  streamwise 


body  surface  location  of  (cp)min  for  each  of  the  bodies,  denoted  earlier 


as  S  ,  is  easily  seen  in  Figure  7  and  shows  that  the  more  blunt  the 
Pc 


min 


body  shape,  the  less  the  extent  of  the  favorable  pressure  gradient  region. 


The  M  curve  can,  to  a  resonably  good  approximation,  be  interpreted 


as  the  6  curve.  The  local  value  of  M  or  6  at  each  point  along  the  body 


surface  can  be  correlated  with  the  slope  of  the  body  shape  obtained  by 
drawing  a  tangent  line  at  the  body  surface.  A  large  slope  gives  a  large 
value  of  6.  This  conjecture  is  based  on  examining  conical  similar  flows 
which  give  constant  values  of  B  for  conical  half  angles  <j> ,  these  values 
of  B  versus  0  are  tabulated  on  page  9  of  this  report.  In  fact,  a  conical 
shape  with  0  large  downstream  of  a  spherical  nose  cap  would  ideally 


represent  a  good  shape  since  B  >  0  over  the  downstream  distance.  However, 


for  Ellipsoidal 


in  order  to  realistically  attach  this  conical  shape  smoothly  to  a 
cylindrical  body  of  prescribed  diameter,  the  conical  surface  distance 
would  be  quite  limited.  Therefore  a  large  pressure  gradient  would  exist 
but  would  not  extend  far  enough  downstream. 

As  previously  discussed,  the  effects  of  body  shape  on  the  growth  of 
the  boundary  layer  thickness  for  a  fixed  does  not  appear  to  be  a 


significant  factor  except  in  the  vicinity  of  (C  )  .  ,  where  an  abrupt 

p  mm 


increase  exists  for  each  of  the  respective  shapes. 

As  emphasized  earlier,  a  large  positive  value  of  8  is  desirable 
downstream  of  the  stagnation  point  for  as  large  a  downstream  distance  as 
possible.  A  tradeoff  clearly  becomes  evident  between  a  large  favorable 
pressure  gradient  (8  large)  and  an  extended  favorable  pressure  gradient 


region  (S 


large) .  A  good  design  of  a  nose  shape  is  one  that  can 


'min 


optimize  these  conditions. 


The  location  of  the  transition  point  St  on  each  body  shape  is  shown 

,6 


in  Figure  7  for  an  =  20  x  10  .  These  results  were  obtained  from  those 
reported  by  Kaups  (MDC  J6530,  April  1974)  in  which  the  Smith-Gamberoni 
prediction  method  was  used.  These  results  confirm  the  arguments  just 
presented.  For  the  elliptic  noses  shown  in  Figure  7,  the  nose  with  an 
aspect  ratio  of  t  =  .50  has  transition  occurring  farthest  downstream. 
Inspection  of  the  8  curves  would  indicate  that  t  ■  .500  and  t  *  .707  have 
the  best  shapes  among  the  four  candidates.  In  summarizing  the  conclusions 


on  the  effect  of  body  shape  on  transition,  one  should  inspect  the 


X  X  ^$3  S 

(C  versus  — ),  (8  versus  r-) ,  and  ( -  versus  — )  curves  for  several  shape 

P  D  D  /R^  D 

noses  and  conclude  that  the  "best  shape"  is  the  one  with  the  following 


attributes : 


J 


#2 


•s 


i)  The  location  of  (C  )  .  should  be  as  far  downstream  of  the 

p  min 


nose  as  possible.  (8  >  0  exists  over  an  extended  region), 


ii)  The  values  of  6  should  be  as  large  as  possible  in  the  upstream 


vicinity  of  S  .  (B  large  positive  number). 

^min  U 

g 

iii)  The  value  of  (— )  should  be  as  small  as  possible. 

U  max 


iv)  The  value  of  should  be  as  small  as  possible  in  the 


upstream  vicinity  of  S 


v)  The  value  of  the  roughness  height  should  be  less  than  20  per 


cent  of  6*  in  the  upstream  vicinity  of  S 


The  conflicting  conditions  (i)  and  (ii)  are  a  direct  consequence  of  the 


fact  that  stability  is  delayed  by  a  favorable  pressure  gradient.  The 


third  conjecture  is  a  result  of  the  fact  that  achieving  (R„  )  * 

us  st  max 

6  t 

( - )  does  not  necessarily  assure  you  that  (S  )  has  been  achieved 

v  max  J  }  t  max 


since  the  local  velocity  U&(s)  is  also  increasing  in  the  favorable  pres¬ 


sure  gradient  region.  Furthermore,  the  boundary  layer  growth  rate 


increases  with  increasing  —  .  Thus,  condition  (iii)  also  introduces  a 


conflict  with  the  criteria  suggested  under  conditions  (i)  and  (ii). 


Therefore,  an  optimum  condition  must  exist  on  choosing  the  shape  with 


rather  unique  pressure  distribution  features. 


Conditions  (iv)  and  (v)  are  based  on  examining  the  results  from  the 


boundary  layer  computations  and  are  of  secondary  importance  in  assessing  a 


good  shape.  Condition  (v)  is  only  a  reminder  that  with  decreasing  6*, 
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desirable  for  delaying  transition,  one  is  more  likely  to  trip  transition 


by  means  of  a  surface  proturberance. 


3.2  EFFECT  OF  HIGH  REYNOLDS  NUMBER  ON  TRANSITION 


In  designing  the  best  possible  nose  shape  for  high  Reynolds  number 
operation  (Rp  >  30  x  10^)  it  is  the  authors'  conjecture  that  the  best  that 


m 
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one  can  hope  for  is  that  transition  will  occur  at  or  slightly  downstream 


of  the  location  of  the  minimum  pressure  (C  )  .  or  that 

p  min 


=  S 

t  P, 


min 


For  these  high  Reynolds  numbers,  the  dimensionless  boundary  layer  rate  of 
growth,  given  by  equation  (A  9 ) 


U 


*d 


increases  as  VR^  and  this  is  illustrated  schematically  in  Figuie  4.  It 

is  clear  that  with  increasing  body  Reynolds  number,  transition  will 

continuously  move  upstream  and  eventually  occur  within  the  favorable 

pressure  gradient  region.  This  is  confirmed  by  the  results  obtained  by 

Kaups  and  illustrated  in  Figure  8  for  elliptical  shape  noses.  At 

Reynolds  number  above  R^  =  30  x  10^,  transition  moves  upstream  of  the 

S  location  for  most  the  bodies  shown  in  Figure  8.  Therefore,  a  large 
min 

value  of  8  upstream  of  the  S  location,  will  slow  down  the 

^min 

upstream  movement  of  the  transition  point  with  increasing  speed  or 
Reynolds  number  R^. 

Furthermore,  the  "best  shape"  for  transition  postponement  at  a  low 
Reynolds  number  might  not  be  the  best  choice  at  higher  Reynolds  number. 
Examining  Figure  8  shows  that  the  elliptical  nose  with  an  aspect  ratio 
of  t  ■  .500  delays  transition  the  farthest  downstream  at  ■  10  x 


10 


while  the  more  blunt  nose  body  with  an  aspect  ratio  of  t  *  .707  delays 
transition  the  farthest  at  a  body  Reynolds  number  of  R^  ■  50  x  10^. 


It  should  be  realized  that  with  an  increasing  R^,  the  thickness  of 
the  boundary  layer  actually  decreases.  This  can  be  deduced  from 
equation  (A  8)  which  gives 


,rd  U 


u 

—  p2 

U  3 


6*(6,0  =  ~ 


8 

_1 

_D _ 

U 
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where 


<5*.(°.--5j -  .568  and  1.2  >  6*(8,0  >  .804 

VT  ~ 


Thus,  the  possibility  of  roughness  tripping  the  laminar  boundary  layer 
increases  significantly  at  high  Reynolds  numbers. 

The  remarks  made  on  the  effects  of  high  Reynolds  number  on 
transition  should  be  interpreted  as  being  somewhat  speculative.  Experi¬ 
mental  data  is  seriously  lacking  in  both  understanding  the  complex 
transition  process  under  these  conditions,  and  confirming  the  validity 
of  the  stability  prediction  method  at  these  high  Reynolds  numbers.  Even 
the  most  basic  types  of  laboratory  experiment  are  not  usually  possible 
in  the  range  of  Reynolds  number  under  consideration.  Secondary  effects 
such  as  surface  roughness  can  play  a  governing  role  on  the  transition 
process  at  these  conditions.  Experimental  tests  on  large  axisymmetric 
bodies  are  required  in  order  to  shed  light  and  give  further  direction 
to  the  factors  governing  transition  under  these  extreme  conditions. 


4.0  SUMMARY  AND  CONCLUSIONS 

A  study  was  made  to  examine  and  gain  insight  into  the  inherent 
features  and  characteristics  of  the  laminar  boundary  layer  on  axisymmetric 


bodies  that  are  used  as  a  basis  for  the  Smith-Gamberoni  method  of  pre¬ 
dictions.  The  ability  of  the  method  to  accurately  predict,  for  a  given 
body  shape  and  prescribed  body  Reynolds  number  R^,  the  location  of 
transition  can  only  be  attested  by  comparing  the  predictions  with  experi¬ 
mentally  obtained  results  such  as  those  obtained  by  DTNSRDC.  However, 
accepting  the  condition  that  this  method  reflects  the  current  state-of-the- 
art  for  prediction  techniques,  an  effort  was  made  to  examine  carefully 
those  potential  flow  and  laminar  boundary  layer  characteristics  of  the 
viscous  flow  region  that  are  used  as  input  for  the  complex  stability 
computations,  and  to  find  a  correlation  between  certain  geometrical  and 
hydrodynamic  characteristics  of  a  given  body  shape,  and  to  ascertain  the 
ability  of  that  shape  in  delaying  transition. 

For  purposes  of  explanation,  the  transition  process  is  envisioned 
as  the  condition  under  which  the  downstream  growth  of  the  boundary  layer, 
denoted  by  R^  ,  is  such  that  it  exceeds  a  critical  value  (R^)  and 
that  the  sophisticated  stability  computation  program  is  essentially  per¬ 
forming  this  computation  as  rigorously  and  accurately  as  possible.  Based 
on  this  interpretation,  the  key  to  delaying  transition  is  to  maintain 
as  thin  a  boundary  layer  as  realistically  possible  and  provide  an 
external  flow  environment  with  a  strong  favorable  pressure  gradient. 

It  is  suggested  that  a  comparative  evaluation  of  various  shapes 
be  made  without  extensive  boundary  layer  computations  or  stability  com¬ 
putations  but  by  simply  comparing  the  shapes  of  the  pressure  gradient 
parameter  B  curves  from  the  stagnation  point  downstream  until  the  C 

^min 

location  is  reached.  It  is  concluded,  with  a  certain  degree  of  con¬ 
jecture,  that  the  best  shape  for  delaying  transition  is  correlated  with 


shape  of  the  8  curve.  It  should  be  emphasized  that  one  can  compute  the 
$  curves  for  each  body  directly  from  the  potential  flow  computations. 

The  question  of  obtaining  the  "optimum  design"  of  a  nose  shape  is  not 
straightforward  but  the  inherent  boundary  layer  characteristics  most 
favorable  in  delaying  transition  are  discussed  and  some  examples  of 
geometrical  shapes  for  which  stability  computations  have  been  performed 
are  presented  to  show  the  attributes  of  a  good  design.  These  examples 
give  some  direction  to  future  work  in  obtaining  optimal  designs. 

An  important  area  for  immediate  assessment  is  the  effect  of 
high  body  Reynolds  number  on  the  transition  process.  This  has  been 
addressed  with  a  certain  degree  of  speculation  since  experimental  data  is 
seriously  lacking  in  either  shedding  light  on  the  transition  process  or 
in  confirming  the  validity  of  the  stability  prediction  methods  at  these 
Reynolds  numbers. 
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For  both  cases  the  following  equations  are  utilized 
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